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A  nonlinear  theory  of  wave  resonances  in  a  two-dimensional  ion  layer  • 
confined  under  the  surface  of  liquid  helium  is  presented.  The  ion  layer  is 
modelled  as  a  two-dimensional  cold  plasma  fluid.  In  addition  to  the  usual 
nonlinearities  present  in  the  continuity  equation  and  the  equation  of  motion, 
the  theory  considers  a  nonlinear  dependence  of  the  mass  of  a  plasma  particle 
on  its  velocity,  as  suggested  by  indirect  experimental  evidence.  Secular 
perturbation  theory  is  used  to  find  the  plasma  response  when  the  damped, 
nonlinear  system  is  driven  externally.  For  typical  experimental  parameters, 
the  mass  nonlinearity  is  found  to  be  the  dominant  nonlinear  effect,  giving 
rise  to  a  backbending  of  the  resonance  curve.  __  _ 
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I.  INTRODUCTION 


Helium  ions  can  be  trapped  just  below  the  surface  of  superfluid  helium. 

These  ions  form  an  almost  ideal  two-dimensional  layer  held  in  position  by 
external  electrostatic  fields.  Waves  can  be  excited  in  the  ion  layer  by  an 
oscillating  electric  potential  applied  to  the  walls  of  the  confining  cell. 

As  the  amplitude  of  the  applied  potential  is  increased,  the  waves  are 
observed  1  to  display  nonlinear  features  such  as  hysteresis,  indicative  of  a 
backbending  of  the  resonance  curve.  It  is  the  purpose  of  the  present  study  to 
provide  an  analytical  description  of  such  phenomena. 

The  linear  properties  of  waves  in  these  systems  can  be  accurately 
described^  by  modelling  the  ions  as  a  cold  two-dimensional  fluid  plasma.  Such 
a  model  is  used  in  this  study  to  explain  the  nonlinear  behavior.  The  first 
source  of  nonlinearity  considered  here  arises  from  the  nonlinear  terms  in  the 
continuity  equation  and  the  equation  of  motion  (the  ponderomotive  effect  term) 
which  describe  the  fluid  motion  in  the  plane  of  the  ion  layer.  The  second 
source  of  nonlinearity  considered  is  a  dependence  of  the  mass  of  the  plasma 
particle  on  its  velocity.  A  brief  discussion  of  what  a  plasma  particle  is  in 
these  systems  is  presented  next  to  motivate  this  proposed  nonlinearity. 

Helium  is  in  the  liquid  state  at  typical  experimental  temperatures 
(<  0.2  K)  and  no  external  pressure.  The  presence  of  an  He+  ion  in  the  liquid 
polarizes  the  surrounding  helium  atoms  and  provides  enough  local  pressure  (due 
to  electrostatic  attraction)  to  freeze  a  sphere  of  about  25  helium  atoms  around 
each  ion.  Since  this  object  moves  in  a  liquid,  its  effective  mass  is 
further  enhanced^  by  half  the  mass  of  the  displaced  fluid.  Thus  the  effective 
mass  of  a  singly  charged  plasma  particle  is  expected  to  be  approximately  150  a.u. 


-3- 


A  fit  of  experimental  data1  with  linear  wave  theory  yields  similar  values. 
However,  it  has  been  experimentally  observed1  that  the  effective  mass  obtained 
by  the  best  fit  between  theory  and  experiments  increases  with  the  temperature 
of  the  liquid  helium  (in  addition  to  depending  on  the  external  confining 
fields  and  the  ion  density).  For  the  temperature  range  in  kef.  1 
(0.1  -0.5  K) ,  the  dependence  is  nearly  linear.  The  microscopic  reason  for 
this  dependence  is  not  yet  fully  understood.  Nevertheless,  one  can 
hypothesize  that  the  increase  in  mass  m  is  associated  with  the  thermal  motion 
of  the  particle  and  that  an  oscillatory  motion  with  a  velocity  v  (caused,  for 
example,  by  a  wave)  is  equivalent  to  an  effective  temperature  Teff  =  m|v|2/2 
(for  two-dimensional  motion  or  twice  this  value  for  one-dimensional  motion). 
This  leads  to  a  dependence  m  =1110(1  +  8|v|2]  where  the  coefficients  mo  and  8 
depend  on  the  external  holding  fields  and  the  static  ion  charge  density.  The 
value  of  6  in  typical  experiments1  is  on  the  order  of  10“  7  c.m“2s2. 

The  experiments  on  ion  layers  have  typically  employed  a  cylindrical 
geometry  and  efforts  are  now  underway  to  use  a  rectangular  geometry .^In  the 
present  work,  the  simpler  Cartesian  model  is  emphasized  for  clarity  of  exposi¬ 
tion,  but  the  corresponding  results  for  the  cylindrical  case  are  also 
described.  Previous  work^ on  the  nonlinear  waves  in  Cartesian  two-dimensional 
plasmas  makes  idealized  assumptions  on  the  equilibrium  density  profile  and 
uses  boundary  conditions  (satisfied  by  the  wave  potential)  which  are  not  well 
justified;  furthermore,  only  the  ponderomotive  nonlinearity  is  considered. 

The  goal  of  the  present  work  is  to  present  a  rigorous  treatment,  _ 
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including  exact  equilibrium  profiles  and  correct  boundary  conditions,  of  the 
nonlinearities  present  in  the  continuity  equation  and  the  equation  of  motion 
as  well  as  the  mass  nonlinearity.  "" 

The  paper  is  organized  as  follows.  The  Cartesian  geometry  of  the  model, 
and  the  basic  cold  plasma  fluid  equations  satisfied  by  the  ion  system  are 
described  in  Sec.  II.  Secular  perturbation  theory  is  used  to  solve  the 
nonlinear  equations.  Static  equilibrium  described  by  the  zeroth  order 
equations  is  obtained  in  Sec.  III.  The  first  order  equations  comprise  an 
eigenvalue  problem  for  the  wave  potential  and  are  discussed  in  Sec.  IV. 
Solutions  to  these  equations  exist  only  for  certain  eigenvalues  of  the 
frequency.  Included  in  the  analysis  of  the  second  order  equations  (Sec.  V> 
are  the  static  ponderoraotive  potential  as  well  as  the  plasma  response  at  the 
second  harmonic.  Analysis  of  the  third  order  plasma  response  at  the 
fundamental  frequency  is  presented  in  Sec.  VI.  The  homogeneous  part  of  the 
equation  is  identical  to  the  first  order  equation  and  therefore  the  existence 
of  a  finite  solution  requires  that  the  inhomogeneous  part  must  be  orthogonal 
to  the  first  order  solution.  This  condition  yields  the  shift  in  the  linear 
resonance  frequency  caused  by  the  various  nonlinear  terms  considered.  It  is 
found  that  the  dominant  contribution  in  a  typical  experimental  ion  layer 
is  the  mass  nonlinearity  which  gives  rise  to  a  backbending  of  the  resonant 
curve.  The  analogous  results  for  cylindrical  geometry  are  presented  in  Sec. 
VII.  Conclusions  are  presented  in  Sec.  VIII. 
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II.  GEOMETRY  AND  BASIC  EQUATIONS 

The  Cartesian  model  geometry  is  shown  in  Fig.  1.  A  two-dimensional 
plasma  strip,  translationally  invariant  in  y  (i.e.,  perpendicular  to  the  plane 
of  the  paper)  and  of  width  2a  in  the  x-direction  is  confined  at  z  =  d  by 
equilibrium  external  potentials  ,  0  and  ^0  applied  to  the  top,  side  and 

bottom  of  a  confinement  cell  of  width  2L  and  height  h.  The  ions  are  held 
just  below  the  liquid  helium  surface  by  the  combined  effect3  0f  the  external 
holding  fields  and  the  dielectric  polarization  of  helium.  The  ion  layer,  which 
is  almost  ideally  two-dimensional,  lies  at  the  liquid  surface  and  forces 
responsible  for  its  vertical  equilibrium  are  not  explicitly  considered  here. 
Furthermore,  since  the  potential  well  in  z  for  the  ions  is  narrow  and  deep, 
and  since  frequencies  considered  here  are  small  compared  to  the  bounce 
frequency  in  this  well,  vertical  motion  of  the  ions  is  ignored.  Finally,  only 
modes  with  wave  vector  component  ky  =  0  are  considered.  These  are  the 
Cartesian  analogs  of  the  azimuthally  symmetric  modes  measured  in  the 
cylindrical  geometry  experiments  of  Ref.  1.  The  cold  fluid  equations  which 
describe  the  two-dimensional  ion  system  have  the  form 

if +  fx  (°'v)  “  o  ’  <l> 


: 


m  r£+w+vfevl  =  -<1  (♦  +  *eo  +  <frewl 

V2$  =  -4jrqa5(z-d)  , 


(2) 

(3) 


where  V2  *  32/3x2  +  32/3z2  .  The  continuity  equation  (1)  relates  the  rate  of 
change  of  the  areal  charge  density  o(x,t)  to  the  velocity  v(x,t)  in  the  plane 
of  the  charges.  Equation  (2)  is  the  equation  of  motion  of  a  particle  of  mass 
m  and  charge  q  under  the  influence  of  $(x,z,t),  the  potential  generated  by  the 
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plasraa  particles,  <£eo(x,z),  the  external  confinement  potential  and  j>ew(x,z,t), 
the  external  potential  applied  to  excite  waves.  As  indicated  in  the  Introduc¬ 
tion,  the  particle  mass  m  is  assumed  to  depend  on  the  velocity  as  m  =  m0  f  1  +  fj  |  v  |  2  ]  . 
Unlike  Eqs.(l)  and  (2)  which  are  defined  only  in  the  plane  z=d  of  the  charges, 
Poisson's  equation  (3)  is  a  three-dimensional  equation  relating  the  plasma 
potential  $  to  a,  subject  to  the  boundary  condition  $  =  0  on  the  walls  z  =  0,h 
and  x  =  ±L  of  the  confining  cell. 

Before  solving  Eqs.  (l)-(3),  it  is  convenient  to  express  them  in  terms  of 
dimensionless  variables  defined  as  follows: 


$  ♦  $ 
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where®  u>nz  =  - 
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with  oo(0)  being  the  equilibrium  density  at  x  =  0.  In  terms  of  the  new 
variables,  Eqs.  ( 1 ) — ( 3)  take  the  form 


I F  +  k<">  "  0 


fl  +  ® I v |  2 1  f  |£  +  vv  +  v|xvl  =  -  (♦  +  *eo  +  few)  z  =  d  ,  (7) 
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72<j>  =  -afi(z-d)  .  (8) 

Secular  perturbation  theory  i.r  used  to  solve  Eqs.  (6)-(8).  The  plasma 
variables  $,  a  and  v  are  expanded  in  a  perturbation  series  with  e  =  5o/oa  as 
the  small  expansion  parameter;  aQ  is  the  equilibrium  density  and  6a  the 
density  perturbation  produced  by  a  wave  driven  by  $ew  =  <j>ej  coscut.  The 
zeroth  order  terms  of  the  expansion  series  represent  the  time-independent 
equilibrium  values  and  the  first  order  terms  represent  linear  waves  with  a 
time  dependence  —  cos  cut.  In  the  absence  of  damping  and  nonlinearities,  the 
plasma  response  has  resonances  at  certain  discrete  values  cuyn  (n  =  1,  2,  .  . 

.)  of  the  frequency  as  shown  in  Sec.  III.  The  inclusion  of  damping  results 
in  a  Lorentzian  resonance  peak  centered  at  tuon,  of  width  v  and  amplitude  4>el/v 
(for  6$).  This  implies  that  for  the  perturbation  expansion  to  be  consistent, 
one  must  have  $el/v<fro  ~  0(e),  where  =  moL2cup2/q .  Including  just  the 
nonlinearities  in  Eqs.  (6)  and  (7)  gives  rise,  in  the  second  order,  to 
time-independent  terms  as  well  as  second  harmonic  terms  ~  cos2u>t  which  provide 
a  third  order  correction  to  the  equation  satisfied  by  terras  proportional  to 
coswt.  This  has  the  effect  of  causing  a  shift  ~  0(e2cuQT1)  in  the  value  of  the 
resonance  frequency  WQn.  When  the  damping  and  the  nonlinearities  are  both 
effective,  the  resonance  curve  has  a  finite  amplitude  and  is  asymmetric  about 
won.  To  mathematically  treat  the  two  effects  together,  it  is  convenient  to 
formally  order  v  ~  0(e2ai0n)  so  that  the  damping  and  the  nonlinearities  first 
appear  in  the  same  order,  namely  the  third;  this  also  implies  that  d>e  1  / d>0  ~ 
0(e3).  The  resulting  approximation  to  the  solution  can,  however,  be  used  for 
all  values  of  v. 


With  the  scales  $ej  ~  0(e3),  the  perturbation  series  tor  a 
and  v  have  the  form 

4>  +  <J>eo  +  4>ew  =  $oo  +  $eo  +  f  1  le  ^  +  c.c.  ]  +  $20 

+  [<j>22e  2l<i)C  +  c.c.l  +  f  ($31  +  $el)e  ^  +  c.c.]  +  (9) 

r  —io)t  i  r  -2iu)t  . 

a00  +  I  <?l  le  +  c.c.  |  +  a 20  +  |  a22e  +  c.c.  | 

+  f °3le  ia,t  +  c.c.l  +  •  •  •  ,  (10) 

c  -iiot  ,  ,  -2ia»t  ,  r  -iu>t  , 

Ivfie  +  c.c.  |  +  |v22e  +  c.c.  |  +  |v3je  +  c.c.  | 

+...,  (11) 

where  vqq  and  V20  are  zero  since  there  is  no  steady  drift  of  particles  in 
the  x-direction.  The  components  on  the  right-hand  sides  of  Eqs.  ( 9 ) — ( 11)  have 
two  subscripts,  the  first  one  referring  to  the  order  of  perturbation  and  the 
second  to  the  harmonic  content.  It  should  also  be  noted  that  the  spatial 

dependence  in  Eq.  (9)  is  on  x  and  z  while  o  and  v  of  Eqs.  (10)  and  (.11)  are 

defined  only  in  the  plane  z=*d  of  the  plasma  layer. 

Equations  ( 9 )— ( 11)  are  substituted  in  Eqs.  ( 6 ) — ( 8 )  and  terms  with  the 
same  time  dependence  and  the  same  order  in  e  equated.  Since  the  modifications 
to  the  resonance  curve  due  to  nonlinearities  is  the  primary  concern  of  the 
present  work,  the  frequency  ui  is  written  a>  =  +  5u>  where  is  the  linear 

resonance  frequency  and  <5u>  =  u>-u>0  ~  0(  v)  ~  0(e2ufc).  The  next  four  sections 
discuss  the  resulting  equations  in  the  first  four  orders  of  perturbation. 
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III.  EQUILIBRIUM 

In  the  zeroth  order,  the  equations  are  time  independent  and  describe  the 
equilibrium  density  profile.  In  addition  to  the  condition  vuo  =  Q,  Eqs.  (b)- 
(8)  in  the  zeroth  order  yield 


f^oo(x,  z  =  d)  +  $eo(x,z  =  d)  j  =  0, 


within  the  plasma,  and 


r32  32  . 

ITT  +  TT'  *oo(x,z)  =  -  a00(x)  6(z-d ) 


3x*  3 z‘ 


subject  to  the  boundary  condition  $00  =  0  at  the  cell  walls.  The  external 
potential  <J>eo  takes  on  the  values  $t0»  0  and  ^0  on  the  top,  side  and  bottom 
of  the  confinement  cell. 

Equations  (12)  and  (13)  can  be  solved  for  the  equilibrium  density  profile 
a00(x).  Using  the  Green's  function,  Eq.  (13)  yields 


(>oo(x,z  =  d)  =  /  dx'cr00(x') 


j  sinh(n  +  -^ ) ird  sinh(n  +  -^l^^h-d)  ^  ^ 

j  - Y  - j -  cos  (n  +  —)irx  cos(n  +  y)irx' 

n  (n  +  y)n  sinh(n  +  -^irh 


Substituting  this  result  in  Eq.  (12)  and  using  a  discrete  grid  for  x  reduces 
the  determination  of  (J00(x)  to  a  matrix  inversion  problem  which  can  be 
uniquely  solved  for  any  choice  of  the  plasma  width  2a  and  the  ratio  <}>t0/$b0* 
The  equilibrium  profile  <J0o(x),  numerically  obtained  for  the  scaled  values 
d  *  0.1,  h  =  0.2,  a  =  0.78,  $tO/$bO  =  ~7  and  for  3  grid  with  256  points  is 
displayed  in  Fig.  2.  The  profile  is  nearly  rectangular  with  the  edge  becoming 
sharper3  with  decreasing  values  of  h. 
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IV.  LINEAR  THEORY 


Only  terms  oscillating  at  the  fundamental  frequency  are  present  in  the 


first  order  equations  arising  from  Eqs.  (6)-(8): 


-i“o°ll  +  (a00  v  1 1)  =  0 


‘,ovll=-j7<j>li(x,z=d) 


V2 d> 1 1  =  -  ai  i  6(z-d ) 


These  equations  can  be  combined  to  give 


-oi02  V2<H1  +  [o0o(x)  '  a'x 1  ]  S(z-d)  =  0  , 


which  on  integrating  across  the  plasma  layer  z  =  d,  yields 


3z  d+  3z  d-  1  dx 


+  Hx  l^aoo<x)  37  4>ll(x>z  =  d>  1  =  0  .  ( 


where  satisfies  the  boundary  condition  i  =  0  at  the  walls  of  the  cell. 


Equation  (19)  is  an  eigenvalue  equation  which  can  be  satisfied  only  for 


certain  discrete  values  of  w0.  It  is  solved  by  using  the  expansion 


cos  (n  +  -t)ttx 


<i>ii(x,z)  -  I  Bn - ; - n  7  TT7 

n  sinh(n  +  yjxd  sinh^n  +  -g- ) n— d ) 


sinh(n  +  i)xd  sinhfn  +  -^)ir(h-z)  ,  z  >  d 


sinh(n  +  yjirz  sinh(n  +  — )x(h-d)  ,  z  <  d  ,  (20) 
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and  the  orthogonality  of  sines  and  cosines  to  yield 

D(o>0)  B  =  0  ,  (21) 

where  the  elements  of  the  matrix  D  are  given  by 

^  sinh(n  +  y)Trh 

^mnC^o)  =  uo2  (n  +  2")  "  JT  ”  K  6mn 

sinh(n  +  J ird  sinh(n  +  — Jir(h-d) 

1  1  ^  1  1 
-  (n  +  -j)  (m  +  -j)  n2  /  dx  aOQ(x)  sinfn  +  ■j)irx  sinfm  +  yjirx  ,  (22) 

and  B  is  a  vector  whose  components  are  Bn  of  Eq.  (20).  The  eigenvalues  u>01, 
u)02»  “03*  •••  for  which  |D|  =  0  and  the  corresponding  eigenfunctions  are 
obtained  numerically.  As  an  example,  the  shape  of  $n(x,z  =  d)  corresponding 
to  the  eigenvalue  ugi  =  0.83  for  the  profile  aoo(x)  is  displayed  in  Fig.  2. 

The  oscillations  seen  in  near  the  wall  are  of  numerical  origin.  The 
x-dependence  of  the  eigenfunctions  on  the  z  =  d  plane  can  be  closely 
approximated3  by  a  single  cosine  function  inside  the  plasma.  The  approximate 
wave  function  cos  3.8x  is  shown  as  a  dotted  curve  in  Fig.  2  for  comparison. 

Since  Eq.  (18)  is  homogeneous,  the  amplitude  of  <j>n  is  undetermined 
within  the  perturbation  scheme.  By  demanding  the  existence  of  a  finite 
solution  to  the  third  order  equation  one  obtains  a  relation  (resonance  curve) 
between  the  amplitude  of  and  the  frequency  w  for  a  given  external  driver 


V.  SECOND  ORDER  EQUATIONS 


The  second  order  equations  contain  the  static  (zero  trequency)  nonlinear 
distortion  as  well  as  the  second  harmonic  response  at  2w  The  zero  trequency 
components  of  Eqs.  (6)-(8)  are 


1 1 v 11* 


+  ffn*vn) 


0 


(23) 


dv  1 1*  dvli  d 

V11  “d 7~  +  vn  “dJ“  =  -  d7  *20(x,z  =  d) 


(24) 


72$20  =  ~  O20  5(z-d)  •  (25) 

Equation  (23)  is  trivially  satisfied  for  and  vn  given  by  Eqs.  (15)— (17). 
Equation  (24)  relates  two  perfect  differentials  and  can  be  integrated  in  x  to 
yield  the  "ponderomotive"  potential 

$20(x,z  =  d)  -  -  | v i i | 2  +  C  ,  (2b) 

where  C  is  an  as  yet  undetermined  constant.  Using  the  Green's  function  tech¬ 
nique  of  Sec.  Ill,  Eqs •  (25)  and  (26)  can  be  reduced  to  a  matrix  equation  and 
numerically  solved  for  020,  given  vn  [from  linear  theory]  and  C.  The 
value  of  C  is  chosen  so  that 

1 

/  dx  o20(x)  -  0  ,  (27) 

-1 

which  implies  that  charge  is  conserved.  The  shape  of  the  nonlinear  static 
density  modification,  o2q(x),  corresponding  to  the  linear  eigenvalue  win  * 
0.83  is  shown  in  Fig.  3.  It  is  observed  that  the  ponderomotive  force  enhances 
the  plasma  density  at  the  center  x=0  and  at  the  edges  x=±a  and  depletes  the 
plasma  elsewhere. 

The  second  harmonic  components  of  Eqs.  (b)-( 8)  are 


V 

•j 

•J 

* 

J 

•J 

•J 
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-  2iu)0  ^22  +  h7  (^oo  v22  +  all  vll>  =  0 


-2iu0  V22  +  vii 


V2^22  =  -  022  5(x-d) 


which  can  be  combined  to  yield 


d 

-  <t)22(x,z  =  d) 


8  3 

(2u0)2  V2(J>22  +  a7  fo00(x)  ^22)  <5(z-d)  = 


3  r  dvn 

“  -57  I  °oo(x)  VI 1  +  2iuio anvil  |  6(z-d) 


Following  Sec.  IV,  Eq .  (31)  can  be  expressed  as  a  matrix  equation 


D^Wq)  e  =  F 


where  the  matrix  D  is  given  by  Eq.  (22)  and  the  vectors  E  and  F  are  defined  by 


*22(x,z  -  d)  -  l  En  cos (n  +  7)1 


"  I7  faooVH  -g“  +  2iuoOHvn]  -  l  Fn  cos  (n  +  y)irx 


The  numerical  solution  of  Eqs.  ( 32)— ( 34 )  for  $22(x>z  =  d)  can  be  readily 
obtained,  once  the  first  order  quantities  an  and  vn  have  been  determined 
from  linear  theory.  For  the  first  few  resonance  frequencies,  2oiQ  is  not  near 
a  zero  of  D(u»)|  and  hence  the  amplitude  of  $22  is  small  (compared  to 
v |J).  rhus,  Eqs.  (2H)  and  (29)  can  be  used  to  obtain  good  approximations  for 
vjj  and  hv  Ignoring  $22  ln  Eq .  (29). 
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VI.  THIRD  ORDER  EQUATIONS 

The  third  order  terms  which  oscillate  at  the  fundamental  frequency 
satisfy  the  equations 


-  i6a)cru  -  10)003!  +  dx  (aoov31  +  <ni*v22  +  02OV  1 1  +  022vll*l  =  0  ’  ('35) 


dvu*  dv22 

-  idow!!  +  wii  -  i^oBlvnl  2vn  -  im^i  +  V22  —  +  vn  = 

~  f$3l(x,z  =  d)  +  *el(x,z  =  d)]  ,  (36) 

V2  <t»31  -  -  03I  «(z-d)  .  (37) 

These  equations  can  be  combined  to  yield 

-  0)o2  v2«J>31  +  1^-  fo00(x)  5(z-d) 

3  3 

“  f_2“o$waii  "  37  fo0  37  Ctel  +  v22vil*^  +  °own 

+  i“o  (®ll*v22  +  O20V11  +  022V11*)  “  1(iloaoB!  vn  |  2vi  1 H  6(z-d)  .  (38) 

The  operator  acting  on  $3!  is  the  same  as  the  linear  Hermitean  operator  of 
Sec.  IV.  This  property  can  be  used  to  find  a  necessary  condition  for  the 
existence  of  a  bounded  solution  $3!  0f  Eq.  (38).  Multiplying  Eq.  (38)  by  <j>n 
and  integrating  over  the  volume  of  the  cell,  it  is  seen  that  the  left-hand  side 
can  be  integrated  by  parts  twice  using  the  boundary  conditions  $n  =  U,  $31  =  0 
at  the  walls  to  obtain  a  volume  integral  of  $32  multiplied  by  the  linear 
operator  acting  on  $n.  But  this  is  zero  from  Eq.  (18).  Therefore,  the  volume 
integral  of  the  right-hand  side  of  Eq.  (38)  multiplied  by  <>32  is  also  zero,  i.e., 


3 

* 

S: 

1 


?: 
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•r, 
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♦>  wvraw 


»v  nu*u  *v* uwv  ■  v*vwv*mw  wfvuii  wwmictiwi 
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-L 


/  dx  4>n(x,z  =  d)  f-  2a)05iooi  i  -  37  [o00  37  |$el(x,z  =  d)  +  v22vn*} 


dx 


+  00W11  +  iu0(ail*v22  +  °20v  1 1  +  a22v  11*) 


“  ia,oa0B|vii|2vnl}  =  0 


(39) 


Writing  =  A  $n  where  is  the  linear  potential  normalized 

such  that  Max.  i|»(x<a)  =  1  [where  ij>  =  aJo_1(d/dx)^n(x,z  =  d)]  and  A  is  the 

complex  amplitude  of  in  this  normalization,  Eq.  (39)  can  be  expressed  as 


“  2  7T  A  “i  +  “2  “  ~  A  ai  +  |A| 2  A  03  =  0  , 


(40) 


where  cq,  a2  and  03  are  real  quantities  defined  as 


1 


oq  =  “o2  /  dx  a00  \Ji2  , 

-1 


(41) 


1 


r  d 

“2  =  «*>0  ]  dx  aOQ  ill  ^7  4>el(x,z  =  d)  , 


(42) 


a3  =  /  dx  {7  a00  t|i  ^7  f ’J'  ^7  Cd>22  -  J  'I'2)! 


-1 


■fix  (*<>♦)  fx  (d>22  -  7  +  2a20’l’2 


dx 


+  I  *2  ^22  "  7  't'2^  "  ♦  77  C°6*)l 


dx 


-  “o2  <*ooB  i))4} 

-  °31  +  <»32  +  <*33  +  034  +  035 


(43) 


with 
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■  $n<x*2  -d> 


(44) 


As  described  in  the  previous  section,  $22  can  be  ignored  (in  comparison  with 
i) i2)  in  Eq.  (43)  for  the  first  few  modes.  Letting  A  =  | A j e i 9  in  Eq .  (40), 
equating  real  and  imaginary  parts  and  eliminating  0  yields 


io  tai  -  «3  t*i*r 


v2cq2 

«i„2 


\  |A| 


2  =  a22 


or 


U» 


"  “o  -  S,  0)0|A|2]"  +  ^  }  |A|  2  = 


2a 


<*22mo  2 

4cq 2 


(45) 


Equation  (45)  gives  the  frequency  response  curve,  | A ]  vs.  w,  where  |a|  is  the 
amplitude  of  for  any  value  of  the  driver  $ei  or  equivalently  a2»  For 

small  values  (i.e.,  «1)  of  the  ’nonlinearity  parameter' 

P  =  (03/2ai)(a22/ai2)(a)03/v3)»  | A |  is  small  and  hence  the  term  (a3u*j/2ai)  |  A|  2 
can  be  ignored  leading  to  the  damped  linear  result,  namely  a  Lorentzian 
response  curve  with  |  A | 2  «  f(ui-u>0)2  +  v2/4]_1.  For  p  >  1,  the  effects  of 
nonlinearity  become  important,  causing  a  frequency  shift  pv  in  the  position  of 
the  response  curve  peak;  the  curve  bends  forwards  or  backwards  depending  on 
whether  p  (or  03)  is  positive  or  negative.  Figure  4  displays  plots  of 
normalized  square  amplitude  (ai2/a22)(  v2/^2)  |  A|  2  vs.  the  frequency 
difference  (a»-ui0)/v  for  the  values  0,  ±1,  ±5  of  p. 

The  numerically  obtained  values  of  04,  02  and  03  tor  the  lowest  four 
modes  (which  are  of  even  parity  in  x)  and  for  $el(x  =  iL)  =  i,  $ei(z  =  0,h)  =  0 
are  presented  in  Table  I.  The  value  of  8  which  causes  the  mass  non¬ 
linearity  <135  to  cancel  the  other  nonlinearities  (031,  032.  033  and  034) 
is  denoted  by  8crit  and  has  the  value  »I0  for  all  the  modes  shown.  Since  the 


% 
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value  of  8  [scaled  as  in  Eq .  (A)]  estimated  for  typical  experiments  isl  o(105), 
the  effect  of  the  mass  nonlinearity  (035)  overwhelms  the  other  nonl ineari t les . 
Since  035  is  negative,  the  frequency  response  curve  bends  towards 
lower  frequencies  with  increasing  amplitude  due  to  a  decrease  in  the  value  of 
of  the  plasma  frequency  as  the  effective  mass  increases  with  the  amplitude 
of  the  wave;  this  leads  to  a  lower  value  of  the  resonance  frequency. 

Good  approximations  to  the  values  of  the  coefficients  aj  given  in  Table  1 

can  be  obtained  by  noting  that  o00(x)  [Fig.  1]  can  be  approximated  by  a 

step-function  density  profile  of  unit  height  and  width  2a.  Also,  one  can  write 

i{i  =  sinKx  as  a  consequence  of  the  single  cosine  approximation  to  $n(x,z  =  d) 

for  |x|<a,  as  illustrated  in  Fig.  2;  the  best  fit  values  of  K  given  in  Table  I 

can  be  approximated  by  Kn  =  mr/a  (corresponding  to  the  resonance  frequency 

Won)*  This  leads  to  ai  3  uio2a.  Approximating  <)>ei(x,z  =  d)  ~  (4i>w/TT)exp[  (x-L)/h] 

if  a  is  not  too  close  to  L,  the  integral  (42)  can  be  performed  yielding 

-1 

“2  -  (-l)n+1  8*w(Knm0/h)f(1T/h)2  +  K^l  exp[  ir(a-L)/h] .  Also,  as  mentioned  in 
Sec.  V,  <t» 2 2  can  be  ignored  in  comparison  with  leading  to  031  3  032  3 
Kn2a/8,  a33  3  -w02(Kna/2)sinh2Knhf sinh2Knd  sinh2Kn(h-d) l"1 ,  034  3  3Kn2a/8 
and  035/8  3  3o)o2a/4.  The  estimate  for  033  makes  use  of  the  approximation 
^20(x,z  =  d)  =  -vi2  +  const  3  -( 1/ 2)cos 2Knx  leading  to 


cos2Knx 

*20(x,z)  3  -  2sinh2Knd  sinh2K„(h-d) 


{ 


sinh2Knd 

sinh2Knz 


sinh2Kn( h-z ) ,  z>d 
sinh2Kn(h-d ) ,  z<d 
,  (4b) 


and  hence  to 


°20(x) 


,3*2() 

>  3z 


3»2() 

3z 


K.n  sinh2Knh 

sinh2Knd  s inh2Kn( h-d )  "nX 


(47) 
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This  value  for  (^q  is  used  in  Eq .  (43)  along  with  ip  ~  sinKnx  to  obtain  the 

estimate  for  0133. 

The  coefficient  033  represents  the  self-modulation  effect  arising  from 
the  static  ponderomotive  density  modification  020  One  can  also  consider  the 
cross-modulation  effect  of  020  (produced  by  the  nth  mode)  on  the  frequency  of 
the  mth  mode.  To  isolate  this  effect,  it  is  convenient  to  use  the  variational 

principle  expression3  for  u)om2  obtained  by  multiplying  Eq .  (18)  by  $11  and 
integrating  over  the  volume  of  the  cell: 

1  d  2 
2  J  dx  0oo(x)  f^-  *llm(x,z  =  d)l 

“om  =  ~  i  2  *  (48) 

/  dz  f  dx  [V^llm(x,2)l 
0  -1 

If  ®oo  is  replaced  by  o00  +  020n  where  O20n  is  a  small  static  density  pertur¬ 
bation  produced  by  the  nth  mode,  then  the  fractional  change  in  u>om2  is 


«(“>om) 


wom 


1 

/  dx  *20n  *m2 
-1 _ 

1 

/  dx  °oo  hi2 

-1 


(49) 


where  ij^  «  (d/dx)$ilm(x,z  =  d).  Approximating  0OO  by  a  rectangular  profile 
of  width  2a  and  using  the  approximations  ^  *  sinKmx  =  sin(mirx/a)  and  Eq . 

(47)  for  02On»  yields  zero  for  the  numerator  on  the  right  hand  side  of  Eq. 
(49)  unless  m  =  n.  Thus  the  ponderomotive  density  modification  due  to  any 
mode  is  expected  to  have  little  effect  on  the  resonance  frequencies  of  other 
eigenmodes.  In  numerical  computations  using  exact  density  profiles  and 
eigenfunctions,  for  example,  a  10%  density  perturbation  produced  by  the  first 


mode  causes  less  than  1%  change  in  the  resonance  frequencies  of  the  next 
three  modes. 
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VII.  CYLINDRICAL  GEOMETRY 

The  Cartesian  results  of  Secs.  II-VI  can  be  extended  in  a  straightforward 
manner  to  the  azimuthally  symmetric  modes  of  relevance  to  the  cylindrical 
geometry  used  in  the  experimental  arrangement  of  Ref.  1.  One  again  obtains 
Eq.  (45)  with  | A |  being  the  amplitude  of  the  linear  mode  and  cq,  (*2  and  03 
given  by 

1 

“1  =  t»02  f  rdr  0oo(O  <J<2  .  (50) 

-1 


r  a 

J^rdr  ooo(r)  t|i  <j»el(r,z  =  d) 


1 

a3  3  -  wo2  8  /  rdr  o0o(r)  ,  (52) 

0 


where  now  the  lengths  (r,  z,  d,  and  h)  are  scaled  in  terms  of  R,  the  radius  of 
the  cell,  the  density  in  terms  of  a00( 0),  the  frequencies  in  terms  of 

Wp  =  [4nq2ooo(0)/moR] 1/2  and  8  in  terms  of  (Rup)-2*  Also  = 

“ o_1  (d/dr)  <frll(r,z  =*  d).  In  writing  Eq .  (52)  it  is  assumed  that  the  contri¬ 
bution  from  the  mass  nonlinearity  is  the  dominant  effect.  To  find  useful 
approximations  of  cq,  a2»  and  013,  it  is  noted  that  in  the  cylindrical  case3 
also,  0Oo(r)  is  nearly  a  rectangular  profile  of  unit  height  and  width  a  and  ij 1 
*  ^l(j,  r/a)  to  a  very  good  approximation  for  0  <  r  <  a.  Also  <£ei(r,z  =  d)  = 
(4$w/ff)  Io(*r/h)/I0(if/h)  if  a  is  not  very  close  to  1.  In  the  expression  for  <p 
and  $el»  J  and  I  are  the  Bessel  function  and  the  modified  Bessel  function, 

respectively  and  j  are  the  zeros  of  Jj.  Using  these  approximations  yields 
In 


(53) 
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a2 


-4«o*w 


In 


1  ( Tra/h ) 


h  (jln/a)z  +  (ir/h)J  I0(ir/h)  J°(jln) 


(54) 


and 


a 

“3  =  -  <^2  s  /  rdr  Jl4(jln  J)  .  (55) 

o 

The  numerical  values  of  alf  a2.  «3  [obtained  from  Eqs.  ( 53)— ( 55) ]  and 

=  a3O22a>o2/2ai30#w2  (which  is  independent  of  (jJq)  are  displayed  in 
Table  II  for  the  first  four  azimuthally  symmetric  modes  for  typical  experi¬ 
mental  parameters^  h  =  0.2  and  a  =  0.895.  The  nonlinear  frequency  shift  is 
pv  as  in  the  Cartesian  case. 


VIII.  CONCLUSIONS 


Using  a  cold  plasma  fLuid  model  for  the  ion  response,  a  nonlinear  theory 
of  wave  resonances  in  a  two-<  imensional  ion  layer  confined  under  the  surface 
of  liquid  helium,  has  been  developed.  Lindstedt-Poincare  theory  of  secular 
perturbation  is  used  to  calculate  the  nonlinear  response  of  the  plasma  when 
damping  and  an  external  driver  are  both  present. 

It  is  found  that  the  usual  nonlinearities  associated  with  the  continuity 
equation  and  the  equation  of  motion  give  rise  to  a  positive  nonlinear  shift  in 
the  value  of  the  resonance  frequency.  However,  the  experimentally  measured 
frequency  shifts  are  negative. 

A  survey  of  the  experiments  suggests  the  possibility  of  another  source 
of  nonlinearity,  namely,  a  quadratic  dependence  of  the  effective  mass  of  a 
plasma  particle  on  its  velocity.  For  instance,  it  has  been  observed1  that  the 
effective  mass  of  an  ion  increases  roughly  linearly  with  the  temperature  of 
liquid  helium,  from  a  value  ~35.2  mne  at  0. IK  to  ~39.6  mHe  at  0.4  K  for  an  ion 
density  of  7xl07  cm-2.  Postulating  that  the  oscillatory  motion  of  an  ion  due 
to  a  wave  corresponds  to  an  increased  effective  temperature  with  AT  =  m|v|2 
(for  one-dimensional  motion),  one  obtains  m  =  m^  [ l  +  g|v|2]  where  mo  =  34  mHe 
and  8  *  6.7xl0~7  cm”2s2.  Scaled  in  terms  of  (Ru>p)~2  where  R  is  the  cell 
radius  and  u>p  =  [4*q2Oo(0)/mR] 1/2,  the  coefficient  8  has  the  estimated 
value  *9x10® .  The  increase  in  the  mass  with  the  amplitude  of  the  wave  leads 
to  a  decrease  in  the  value  of  oip  and  to  a  decrease  in  the  value  of  the 
resonance  frequencies.  The  large  value  of  the  scaled  B  insures  that  the 
negative  shift  of  the  resonance  frequency  due  to  mass  nonlinearity  overwhelms 
the  contributions  from  the  nonlinear  terms  in  the  continuity  equation  and  the 
equation  of  motion.  Using  Table  II  and  the  valuel  v  =  2.74xl0~2  u>p,  one 
obtains,  for  the  lowest  azimuthally  symmetric  mode,  w2^  “L2  (I  -  c  #y2) 


where  $w  is  measured  in  volts  and  Z  2  1100  V-2.  This  value  of  z  is  to  be 
compared  against  the  experimental  value  of  2o0  V~2.  There  is  also  some 
experimental  evidence^  that  v  also  increases  nonlinearly  with  amplitude;  this 
would  lower  the  theoretical  prediction  for  5.  However,  the  important  point  is 
that  our  conjecture  of  a  mass  nonlinearity  gives  the  correct  sign  of  the 
effect  and  yields  results  close  to  the  experimentally  observed  values,  in 
spite  of  various  uncertainties  regarding  indirect  measurements. 

Since  at  the  present  time  there  does  not  exist  a  first-principles  theory 
capable  of  predicting  the  value  of  B,  and  since  it  is  also  possible  that  its 
value  may  vary  from  one  experimental  set-up  to  another,  in  this  study  we  have 
introduced  the  concept  of  a  critical  8.  Its  numerical  value  is  illustrated  in 
Table  I.  Physically,  8Crit  defines  the  regime  beyond  which  the  usual  fluid 
nonlinearities  become  less  important  than  the  mass  nonlinearity.  This  may  be 
a  useful  tool  in  assessing  the  operational  regime  of  future  experiments. 
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FIGURE  CAPTIONS 

Fig.  1  Cartesian  model  geometry. 

Fig.  2  Scaled  equilibrium  density  profile  a00(x)  for  the  scaled  (to  L) 

values  d  =  0.1,  h  =  0.2,  a  =  0.78,  4>to/<J>bo  =  “7  and  the  eigen¬ 
function  $n(x,z  =  d)  (solid  curve)  corresponding  to  the  lowest 
eigenvalue  togi  =  0.83.  The  dots  are  a  plot  of  cos3.8x. 

Fig.  3  Nonlinear  static  density  modification  02q(x)  (produced  by 

of  Fig.  2)  scaled  in  terms  of  100  [Max  v^(x  <  a)]2* 

Fig.  4  Frequency  response  curve  in  which  the  square  of  the  amplitude 

| A |  is  plotted  vs.  the  driver  frequency  for  the  values  0,  ±1,  ±5 
of  the  'nonlinear'  parameter  p  =  (a3/2cti)(a22/ctl2)(|Jl>o3/'>3)* 


i 
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TABLE  CAPTIONS 

TABLE  I.  Scaled  numerical  values  of  uIq,  ft,  oq,  02,  0131,  032.  333,  034 

and  035/8  for  Cartesian  geometry  and  the  scaled  (  to  L)  values  d  = 
0.1,  h  =  0.2,  a  =  0.78  and  $to/$bo  =  -7.  Bcrit  is  the  value  of  8 
for  which  the  frequency  shift  produced  by  the  mass  nonlinearity 
cancels  the  frequency  shift  produced  by  fluid  nonlinearities. 

TABLE  II.  Scaled  numerical  values  of  04,  02,  33  and  the  nonlinearity 

parameter  p  obtained  from  the  approximate  Eqs.  C 53)— ( 55)  for  the 
scaled  (to  L)  values  d  =  0.1,  h  =  0.2,  a  =  0.895  corresponding  to 
typical  experiments  in  cylindrical  geometry. 
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